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Abstract 

<N 

We verify functional a posteriori error estimates of numerical solutions of a contact problem in ID. The 
contact problem is described by a quadratic minimization problem with a bound constrain and it can be 
equivalently formulated as a variational inequality. The simplification into ID allows for the construction 
•^j- of a benchmark for which an exact solution can be derived. Numerical solution obtained by the finite 

element method can be compared with the exact solution and their difference is further estimated in 
the energetic norm by a functional majorant proposed by Repin et al. Numerical tests provide sharp 
estimation of the approximation error. 

1 Formulation of contact problem in ID 



We will deal with a state problem motivated by the following linear boundary value problem: Find a function 
u such that 



^ -u" = f in (0,1) (1) 

m u(0) = u(l) = 0, (2) 

(N 

where / is a continuous right-hand side, / < 0. Equation ([I]) describe, for instance, the elastic deformation 
of a string. Here, / stands for an external force density. If an obstacle (ft is added, see Figure 1, then a 
contact problem for elastic deformation of string arises. This problem is typically described by a variational 
fT) inequality equivalent to the following minimization problem: 

Problem 1 (minimization problem). Find u 6 K satisfying 



l r i 
') 2 < 

JO 



J(v) := i / (t'rd.r - / frd.r [ I) 



and the admissible set is defined as 

K := {v G V : v[x) > <f>(x) a.e. in (0, 1)}, 

where / £ L 2 (0, 1), Vq denote the standard Sobolev space H^(0, 1) and <p £ H 1 ^, 1) such that (/> & V and 
4>{x) < a.e. in (0, 1). 
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Figure 1: Illustration of the contact problem. 

Problem [I] is a quadratic minimization problem with a convex constrain and the existence of its minimizer 
is guaranteed by the Lions- Stampacchia Theorem [3J. It is also equivalent to the variational inequality: 
Find u G K such that 



/ u'[v — u)'dx > / f(v — u)dx for all v G K. 
Jo Jo 



(5) 



The convex constrain v G K can be transformed into a linear term containing a new (Lagrange) variable in 

Problem 2 (Perturbed problem). Let W := {v + tcfi : v G Vq and t G M} C -ff 1 (0, 1), where M. stands for set 
of real number. For given 



H £ A := {fx G W* : (jm, v - <j>) > for all v G K} 



find M„ G Vq such that 



J^(u M ) = inf Jf,(v), 

vEVo 



(6) 
(7) 

(8) 

0) 

(10) 
(11) 

Proof. Obviously, W is the least linear subspace of H l (0, 1) such that definition ^ of perturbed functional 
Jp, is correct. Let w G W is arbitrary. We decompose 



where W* denote the dual space of W and the perturbed functional is defined as 

J M (u) := J(v) - (n,v-(f)). 
Problems [T] and [2] are related and it obviously holds 

Jfj,{Ufi) < J(u) for all n £ A. 
Lemma 1 (Existence of optimal multiplier) . There exists A G A such that 

U\ = U 

and 

J\{u) = J(u). 



w = v + t(f), 



:i2) 



where v G Vo and t G K and this decomposition can be shown to be unique. Now, we define a functional A 
as follows: ^ 

(A, w) := j u'v'dx — / fvdx + t 
Jo Jo 



1 ri 
(u') 2 dj; — / /ndx 
o •/ o 



(13) 



We assert that the functional defined by (13) has required properties (10) and (11). Apparently, A is a linear 



functional on W. The functional A is also continuous. It is a consequence of continuity of decomposition 



(12), which can be proved as follows. Let w G W is arbitrary and w n — > w in W, where w n G W. With 
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respect of ( 12 ), we can write w n = v n + t n cp and w = v + tcft, where v n , v G Vq and t n ,t G R. As a consequence 
of unique orthogonal decomposition of <p we infer 

\\w n - H| H i( 0j i) ^ K - 1 \\\4> ± IIhi(o,i)' 

where (p 1 - is the component of <j> orthogonal to subspace Vq. Moreover, it follows from the triangle inequality 
that 

11% - f|| #1(0,1) ^ \\ w n ~ W\\m(0,l) + K ~ A 11011^1(0,1)- 

Consequently, v n — > v in W and t n — > t. Now, if we restrict the space W to the origin Vq, we obtain 



u'v'dx 



fvdx for all to G Vq, 



(14) 



which is equivalent to the property (10), i.e., inf J\(w) = J\(u). Furthermore, if we take w = u — (f>, it 
follows from (13) that 



(A, u — 4>) = (A, u) — (A, ( 







and consequently the property (11) is fulfilled. Finally, we verify the condition of nonnegativity from 
definition ([6]). Let v G K is arbitrary. It follows from ^ that 

(\,v-4>) = (X,v-u) > 0. 



□ 

Remark 1 (Existence of optimal multiplier with a nonpositive obstacle <j) G Vo). If we would deal with an 
obstacle G Vo, the existence of optimal multiplier could be proved as follows: Once again, the relation (14) 
defines a linear continuous functional A in Vq such that u minimizes the perturbed functional J\ defined by 
@. Since <fi £ K, we can choose v = 4> as we U as v = 2u — (p and apply the inequality ^ , which implies that 
(A, u — 4>) = 0. Subsequently, it follows from ^ that J\(u) = J(u). The condition of nonnegativity from 
definition ^ is also fulfilled, it is sufficient to apply the inequality §5§ to a function u + v, where v G Vq and 
v > a.e. in (0,1). 

Remark 2 (Representation of (13) by a function A G L 2 (0, 1)). If u G -ffo(0, 1)0 -f^ 2 (0, 1), then integration 
by parts yields 



(X,w) 



u"vdx 



fvdx + 1 



u'udx 



fudx 



Xvdx + t Andx, 



where 

X = -(u" + f). (15) 
In this case, the condition (f> G H 2 (Q, 1) have to be fulfilled. 

Remark 3 (Saddle point problem interpretation) . The pair (A, u) is apparently a saddle point of the perturbed 
functional satisfying 

J\(u) = sup inf Jfj,(v) = inf sup J^v) = J{u). 
For more details, you can see e.g. [U p. 107-109]. 



2 Benchmark with known analytical solution 

We derive an exact solution of Problem [T] assuming negative constant functions / and 4> ■ The resulting 
solution is displayed in Figure [2] for the case of active obstacle. A mechanical intuition suggests that for 
small values (considered in absolute value) of acting force /, there will be no contact with the obstacle and 
there will be a contact on a subset of interval (0, 1) located symmetrically around the value x = 1/2 for 
higher values of /. Note that the solution of Problem [T] with inactive obstacle is equivalent to the solution 
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Figure 2: Benchmark setup: Constant forces / pressing continuum against a constant lower obstacle 
exact displacement u (left) and construction of exact displacement u in detail (right). 



of Q - Q and it reads 



u{x) = —(a; — x 2 ). 



(16) 



The minimal value of u on interval (0,1) is attained at x = 1/2 and the inactive obstacle condition u(l/2) > 
is satisfied for 

l/l< 



Then, the corresponding energy reads 

The obstacle is active if 

and the solution has the following form 



J{u) := 
1/1 > 



24' 



(17) 

(18) 
(19) 



l x 2 + <t>+t[\-rf 



u(x) = < 



if x G [0, 



if x £[\-r,\+r] 



(ar— 1) if K€(|+r,l] 



for unknown parameter r G [0, i]. The parameter r determines the active contact set — r, ^ + r] and its 
value can be determined from the minimum of energy 



J(u) 



■ + i(i-r) 2 ] a 



Z/ 1 



2[0 + f (3 " 01/(2 " + 3 ( 2 " r) " " 2M 



over all value of r G [0, i]. The minimal energy 



is achieved for the argument 



1 

r = - 
2 




(20) 



(21) 



Therefore, the solution of the problem with the active obstacle reads 

if x G 
if x G 



ii(a?) = < 



■£x 2 - ^Wf 



X 



'M 1 



(22) 



-{(x-l) 2 + v / 2W(x-l) if 16(1-^,1 



Figure [3] provides few numerical approximations of u, see Section [5] for details. The first-order derivative 



4 



lower obstacle 
solution for f=-6 
solution for f=-8 
solution for f— 10 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



Figure 3: Solutions for problems with loadings / G {—6, —8, —10} and = —1. 



u'(x) 



-fx-y/m 



if x g 

if x G 



0: 



(23) 



-f{x-l) + y/2$J if a?e (1 

is continuous everywhere. It is not difficult to show that 

-/ if x G (0, M 



u"(x) 



I if XG (1_ /Ml 



(24) 



is the second-order weak derivative of (22). Thus, it follows from (15) that the optimal multiplier for our 
benchmark problem reads 

' if x G (O, J?f 



X(x) 



so that it is a piecewise constant function. 



-/ if x G 



if x G I 1 



1 

/ ' 1 



(25) 



3 Functional a posteriori error estimate in ID - general part 

We are interested in analysis and numerical properties of the a posteriori error estimate in the energetic 
norm 



v E ■-- 



(v') 2 dx 



y 




Figure 4: Contact zone of the discrete solution v covers the contact zone of the exact solution. 



This section is based on results of P. Neittaanmaki and S. Repin [5]. We apply the general framework to 
one dimensional case, but try to achieve exact integration rules wherever it is possible. As well as in [5], it 
is simple to see that 



J(v) — J(u) = — I [(v — u)'] 2 dx + I u'(v — u)'dx — f f(v — u)dQ 
2 Jo Jo Jo 



and ([5]) implies the energy estimate 

1 



v — U \\E — J( v ) ~ J( u ) f° r an v E K. 



(26) 



(27) 



Estimate (27) can only be tested for problems with known exact solution u S K. One problem with known 



exact solution will be given in Section [2j 



Remark 4 (Sharpness of estimate (27)). It is clear from (26), the estimate (27) turns into equality if 



(A, v — u) = I u'(v — u)'dx — / f(v — u)dx = for all v € K. 
Jo Jo 



(28) 



This situation always occurs if A = 0. Then, (14) implies that u is a solution of Problem [T] in the whole 
space Vq. This corresponds to a linear problem without any obstacle. However, the estimate (|27|) can turn 



into equality also for the active obstacle. Indeed, (28) rewrites as 

(X,v- u) = -f 



1- ^ 
V f 



(v — u)dx = 0, 



if the contact zone of an approximate solution v G K includes whole contact zone 
solution u, see Figure |4j 



(29) 



of exact 



Further, we substitute (|9|) into (27) to obtain the estimate 



J(v) - J (it) < J(v) - J^Up) for all \x 6 A. 



(30) 



In practical computations, £ A will be approximated by u^h G Vo,h from some finite dimensional subspace 
Vo,h C Vq. Therefore, it holds 



and Ju(uu,h) can not replace J^{u^) in (30) so that the inequality holds. Therefore, estimate (30) is 



practically not interesting. To avoid this difficulty, we establish the following problem: 
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Problem 3 (Dual perturbed problem). Find r* G Q*^ C L 2 (0, 1) such that 

j;{r;)= sup j;( q *), 



where 



and 



JZ(q*) = -~l (<?*) 2 d* + (/v 



Q}^ := I q* G L 2 (0, 1) : = gVdx - J fvdx for all u G Vb| 



sup J*(<f ) = J*(«' M ) = J M («,«). 



Lemma 2. It holds 

Proof. As a consequence of (p7j) , it holds that G Q*f fl - Let u; G is arbitrary. Since 



/a* 



1 



(31) 

(32) 
(33) 



and J Q u'^w — u'^)dx = in consequence of u'^ w G Q*f^, we deduce that J*^{u'^) is supremum of dual 
perturbed functional J* Finally, it is not difficult to verify that Ju(u') = J^u^). 

□ 



Corollary 1 . If we choose 
it holds 



M = A, 



J A (u) = inf J x (v) = sup J* x (q*) = J* x (u') = J(u). 

^0 g *gQ* A 



It follows from Lemma 2l we can replace the inequality (30) by 



J(v) - J(u) < J(v) - sup j;(q*) < J(v) - j;(q*), 



(34) 



where q* G Q*j^ is arbitrary. It is easy to see that we can rewrite ( 34 ) to form 



J{v) - J{u) < 



\(v'?-fv 



dx + J' \(r*) 2 dx +IJ 1 [(q*) 2 ~ (r*) 2 ] dx - (/i, 



(35) 



where r* G L 2 (0, 1) is arbitrary. 

From now, we only consider more special case: 

H e A : = [y e L 2 (0, 1) : y, > a.e. in (0, 1)} . 



Then, 



l 

(/x, v) = I yvdx. 
'o 



The application of estimate (35) is possible, the only limitation is to satisfy the constrain q* G Q*f fM - In the 
following, we recall the technique of [5] in order to replace the constrain by a special penalty term. We 
obtain 



J(v) - J{u) < 



V+ 1 f \v' -T*) 2 dx+^[l + -\ I (q*-T*) 2 dx+ I lx(v-<t>)d X , 



1 



P) Jo 



where j3 > 0. 
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Since q* G Q*^ is arbitrary, we can take the error estimate 

^ J\v' ~ r*) 2 dx + \ (\ + 1) [d 2 (r*, Q/ M )] + jf* /i(t; - </>)dz, 



J(v) - J{u) < 

where 

d 2 (r*,Q* f „) := 

Jo 



d\r*,Q%):= inf [\q*-T*fdx. 



If (r*)' G L 2 (0, 1), i.e. t* G fi^O, 1), we can infer, analogously as in [5j, that 

d 2 (.T*,Q%)<C 2 \\(T*y + f + v\\h {n) , 
where a constant C > originates from the Friedrichs inequality 



f u 2 dx < C 2 [ (u') 2 dx Vu G £#(0, 1). 
Jo Jo 



Since, in our problem, we can put C = 1, we obtain ultimately the following error estimate: 

J(v) - J{u) < I±£ j\v> - r*) 2 dx + \ (\ + ||(r*)' + / + M ||| 2(n) + jf 1 M (« - 0)dx, (36) 



where the right-hand side of (36) is so called majorant M. = M. (v, /, (ft; (3, /i, r*). 



Remark 5 (Optimal majorant parameters). Let fj, = X, t* = u' and /3 — >■ 0. Then, it is not difficult to 



show that the inequality in (36) changes to equality, i.e. the majorant on right-hand side of (36) defines the 



difference of energies J(v) — J(u) exactly. 

4 Majorant minimization 

For given solution approximation v , loading / and the obstacle (ft, the majorant M represents a convex 
functional in each of unknown variables /3, /i ,r*. Our goal is to find such variables /3 op t, ^opt and r* t that 
minimize the majorant Ai. 

Problem 4 (Majorant minimization problem). Let v G K, f G L 2 (0,1), (/> < be given. Find optimal 
Popt > 0, /x opt G A and r* pt G -ff 1 (0, 1) such that 

M(v, f, (ft; /3 opt , fJ- op t,T* t ) = min M(v, f,(ft; /3, fj,,T*). 

l3,fl,T* 

Algorithm 1 (Majorant minimization algorithm). Let k = and and let > and /ij. 6 A be given. Then: 
(a) find t% +1 G -£^(0, 1) such that 



M(v,f,(ft;/3 k ,n k ,T k+1 ) = min M(v, f, (ft; f3 k , (j, k , r*), 

T*ei/ 1 (o,i) 



(b) find /Ufc + i G A such that 



M(v,f,(ft;P k ,fi k+ll Tl +1 ) = minAl(u,/,0;/3 fc)) u,r| +1 ), 



(c) find /3fc + i > such that 



, /, (ft; P k +\ , ^fc+i , T k+i ) = M (v, f,(ft; /3, n k+x , T k+1 ) , 



(d) set k = k + 1 are repeat (a)-(c) until convergence. 
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The minimization in step (c) leads to the explicit relation 

Pk+i = 



T k+l)' + / + M*+i||l2(o,i) 



\v -r fe+1 || L 2 (0jl) 



(37) 



The minimization in step (a) is equivalent to the following variational equation : Find t^ +1 £ -ff 1 (0, 1) such 
that 



1 + p k ) v'wdx - ^1 + ^\ J (/ + m,)v/dx for all w G ^(0, 1) 



(38) 



In later numerical experiments, we search a solution of (38) in a finite-element space Vh C ff 1 (0, 1). By 
introducing an interval partition with nodes 



= X\ < X2 < ■ ■ ■ < X n = 1, 



(39) 



we define Vh as the space of functions continuous on [0, 1] and linear on the interval [xj, Xj+{\, j = 1, . . . , n— 1. 
In usual nodal basis ipj, j = 1 . . . n, the relation (38) rewrites as a linear system of equations in matrix form 



;i + p k )M + ( i + — ) a 



y={l + Pk)b - ( 1 + — |r. 



(40) 



for unknown coefficients y = (yi, . . . ,y n ) G i? n representing r^ +1 = Y^j=iVi' l l ) r Here, M = (m^) is a mass 
matrix, ^4 = (ay) is a stiffness matrix, 



rriij = j ipiipjdx, dij = / tp'iip'jdx 
o Jo 



and & and c are n- dimensional vectors 

h = 

for i,j = l... n. 



1 Z" 1 

v'tpidx, Q = / (/ + ^kWidx 
o jo 



The minimization in step (b) is equivalent to the variational inequality: Find /ifc+i £ A such that 



l + J k ) [Mk+l+(-A- + l)' + /] + 



(io — jUjt + i)dx > for all mi£ A. 



In numerical experiments, we look for fik+i € A in the space of piecewise constant functions defined on (39). 



Under the assumption of piecewise linear functions 4> and v and piecewise constant functions / and (r*) on 



( 39 ) with given values 



(Xj),<f>(x j+1 ), v(Xj),v(xj+l), f(x j+ l), (T% +1 )'(x j+ l), 



the exact minimizer reads 



Vk+l(x j+ i) = -(T% +1 )'(x j+ i) - f{x j+ i) 



v(Xj) + v{Xj+l) - (j)(Xj) - (/>(Xj +1 ) 



2 ^ + m 



(41) 



where (-) + = max{0, •}. 
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5 Numerical experiments 



Numerical approximation v of the exact solution u from (22) is searched as a piecewise linear function 



defined on the partition (39) and constructed by using the Uzawa algorithm. 
Algorithm 2 (Uzawa algorithm). 

1. Set the initial Lagrange multiplier p,Q = 0. 

2. Start of the loop: for /c = 1,2,... do until convergence: 

3. Find v k G Hq(0, 1) such that J^Vk) — > min. 

4. Set Hk = (fJik-i + p(vk ~ 4>)) + - 

5. End of the loop. 

6. Output v = v k and p = p k - 

The minimization in step 3. of Algorithm [2] is equivalent to the following variational equation: Find 
v k G Hq(0, 1) such that 

/ v' k w'dx = I (/ + pk)wdx for all w E Hq(Q, 1). 
Jo Jo 



By using the partition (39) and corresponding nodal basis, we can represent the k-th iteration v k by linear 
combination of the basis functions 1^2, ■ ■ ■ , VVi-l? i-e. v k = S?=2 %" i Ar Then, the unknown coefficients 



2/2> • ■ • j Vn-1 can De obtained by solving of a linear system of equations analogous to (40). 



The convergence of Algorithm [2] depends on the choice of the scalar parameter p and it can be shown, 
see e.g. [?], that is alway converges for p 6 (0, p\) for some p\ > 0. Here, we consider the case / = —14, 
choose p = 10 and run Algorithm [2] with 10000 iterations. The first iteration (vi,p±), the second iteration 
(v2,p-2) and the final (the 10000th) iteration (v,p) are displayed in Figure [5j 

Assuming the uniform mesh with 641 nodes (which correspond to 6 uniform refinements of an initial 
uniform mesh with 10 elements), we consider approximate solutions (v,p) obtained by Algorithm [2] with 
100, 1000, and 10000 iterations and we verify estimates (27) and (36). The results are reported in Tables [TJ 



[2j[3j For cases with an inactive obstacle (for loadings / G {—5, . . . , —8}), the energy estimate (27) is always 
the sharpest possible (see values 1.00 in the fourth column of each Table), which is known from Remark |4j 
Such cases are described by purely quadratic minimization problems, whose discrete solutions are obtained 
already after one iteration of the Uzawa Algorithm. In other cases with an active obstacle (for loadings 
/ G {—9, . . . , —18}), we a priori choose the numbers of iterations in the Uzawa Algorithm in a way, that: 

100 iterations provide inaccurate approximations v on the given uniform mesh of exact solutions u, 

1000 iterations provide more accurate approximations v on the given uniform mesh of exact solutions u, 

10000 iterations provide very accurate approximations v on the given uniform mesh of exact solutions u. 

The majorant minimization algorithm is run with a high number of iterations (10000 in all experiments) 



in order to achieve the sharpest possible estimate (36) (this corresponds to values around 1.00 in the last 



column of each Table). We note that the functional majorant provide a sharp estimate of the difference of 
discrete and exact energies, no matter what quality of the approximation v above we consider. 

Remark 6 (Update of /3). The experiments showed that the update of j5 in the step c) of Algorithm [Tj should 
not be called in every iteration. It turns out useful to call steps a) and b) repeatedly and run step c) only 
after variables r* and p stabilize. We updated (3 during the 5000th and the final 10000th iterations. 
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The iterations of the majorant minimization algorithm are displayed in Figure [6j We consider initial 
setup Po = 1, fj,o = and the inaccurate approximation v provided by 100 iterations of the Uzawa algorithm 
(taken from above). The majorant minimization algorithm with 10000 iterations provides a good quality of 
/x as well as of t* (see Remark [5]) . 

Here, we are not interested in a fast and efficient majorant optimization, new strategies combined with 
effective minimization known for linear elliptic problems |6] will be developed in forthcoming papers. 
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Figure 5: The first, the second and the final (the 10000th) iterations of the Uzawa algorithm for / = —14. 
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Figure 6: The first, the twentieth and the final (the 10000th) iterations of the majorant minimization 
algorithm for / = —14. in 
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6 44e-004 


2.70 


6 51e-004 
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1.02e-004 


9.91e-004 


3.11 


9.98e-004 


1.00 


-17 


1.13e-004 


1.26e-003 


3.34 


1.27c-003 


1.00 


-18 


1.24e-004 


1.54e-003 


3.53 


1.55e-003 


1.00 



Table 1: Verification of majorant and energy estimates for problems with various / computed on an uniform 
mesh with 641 nodes. Discrete solutions v is computed by 100 iterations of the Uzawa algorithm. 



/ 


\h-u\\l 


J{v) - J{u) 


/J(v)-J(u) 

V Ilk-"!!! 


M(v,...) 


/ M(v,...) 
V J(v)-J{u) 


-5 


2.54e-006 


2.54e-006 


1.00 


2.55e-006 


1.00 


-6 


3.66e-006 


3.66e-006 


1.00 


3.67e-006 


1.00 


-7 


4.98e-006 


4.98e-006 


1.00 


4.99e-006 


1.00 


-8 


6.51e-006 


6.51e-006 


1.00 


6.52e-006 


1.00 


-9 


9.04e-006 


9.49e-006 


1.02 


9.83e-006 


1.02 


-10 


1.00e-005 


2.36e-005 


1.53 


2.39e-005 


1.01 


-11 


1.19e-005 


2.72e-005 


1.51 


2.76e-005 


1.01 


-12 


1.33e-005 


2.90e-005 


1.48 


2.94e-005 


1.01 


-13 


1.53e-005 


3.96e-005 


1.61 


4.01e-005 


1.01 


-14 


1.73e-005 


4.79e-005 


1.66 


4.85e-005 


1.01 


-15 


1.88e-005 


5.61e-005 


1.73 


5.69e-005 


1.01 


-16 


2.07e-005 


5.51e-005 


1.63 


5.59e-005 


1.01 


-17 


2.24e-005 


6.01e-005 


1.64 


6.10e-005 


1.01 


-18 


2.48e-005 


6.94e-005 


1.67 


7.03c-005 


1.01 



Table 2: Verification of majorant and energy estimates for problems with various / computed on an uniform 
mesh with 641 nodes. Discrete solutions v is computed by 1000 iterations of the Uzawa algorithm. 



/ 




JO) - J{u) 


/J(v)-J(u) 


M(v,...) 


/ M(v,...) 
V J(v)-J(u) 


-5 


2.54e-006 


2.54e-006 


1.00 


2.54e-006 


1.00 


-6 


3.66e-006 


3.66e-006 


1.00 


3.66e-006 


1.00 


-7 


4.98e-006 


4.98e-006 


1.00 


4.98e-006 


1.00 


-8 


6.51e-006 


6.51e-006 


1.00 


6.51e-006 


1.00 


-9 


7.78e-006 


8.05e-006 


1.02 


8.10e-006 


1.00 


-10 


9.11e-006 


9.79e-006 


1.04 


9.88e-006 


1.00 


-11 


1.05e-005 


l.lOc-005 


1.02 


l.llc-005 


1.00 


-12 


1.20e-005 


1.29e-005 


1.04 


1.30e-005 


1.00 


-13 


1.35e-005 


1.42e-005 


1.03 


1.44e-005 


1.00 


-14 


1.51e-005 


1.60e-005 


1.03 


1.61e-005 


1.00 


-15 


1.68e-005 


1.78e-005 


1.03 


1.79c-005 


1.00 


-16 


1.84e-005 


2.01e-005 


1.04 


2.03e-005 


1.01 


-17 


2.02e-005 


2.16e-005 


1.03 


2.18e-005 


1.00 


-18 


2.20e-005 


2.39e-005 


1.04 


2.42c-005 


1.01 



Table 3: Verification of majorant and energy estimates for problems with various / computed on an uniform 
mesh with 641 nodes. Discrete solutions v is computed by 10000 iterations of the Uzawa algorithm. 
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